Fluctuations of Fluctuation— Induced "Casimir" Forces 
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The force experienced by objects embedded in a correlated medium undergoing thermal 
fluctuations — the so-called fluctuation-induced force — is actually itself a fluctuating quantity. We 
compute the corresponding probability distribution and show that it is a Gaussian centered on the 
well-known Casimir force, with a non-universal standard deviation that can be typically as large as 
the mean force itself. The relevance of these results to the experimental measurement of fluctuation- 
induced forces is discussed, as well as the influence of the finite temporal resolution of the measuring 
apparatus. 
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In 1948, Casimir predicted that two uncharged con- 
ducting plates facing each other in vacuum are subject 
to a long-range universal attraction Q. The latter is 
due to the modification of the quantum electromagnetic 
fluctuation spectrum due to the boundary conditions im- 
posed by the conducting plates. In the last fifteen years, 
the concept of Casimir force has been extended to ther- 
mally excited elastic forces between objects embedded in 
a medium with scale- free fluctuations altered by the pres- 
ence of the objects [Q. Soft matter examples comprise 
inclusions in complex fluids undergoing thermal fluctu- 
ations Jpj, e.g. fluid membranes || [| and critical mix- 
tures [pL or interfaces bounding complex fluids, e.g., liq- 
uid crystals || Q and superfluids ||, ^| . One of the main 
features of these Casimir interactions is their universal- 
ity. In a given geometry and at a given temperature, the 
fluctuation-induced forces depend only on the universal- 
ity class of the fluctuating medium and on the nature of 
the imposed boundary conditions (the material's elastic 
constants and the coupling strength of the boundaries 
are irrelevant). 

While the development of micro- and nanoscale experi- 
ments has allowed precise direct verifications of the origi- 
nal vacuum's Casimir force [llj , experimental charac- 
terization of thermal fluctuation-induced forces in soft- 
matter has been scarce [|l3| . The main difficulty with 
measuring these forces is usually imputed to the presence 
of stronger background elastic and van der Waals forces. 
Another important point which has been overlooked, is 
that the fluctuation-induced force is itself a fluctuating 
quantity: what is commonly referred to as the "Casimir 
force" is only its ensemble average. The fluctuating na- 
ture of this force should be taken into proper account 
before their complete characterization can be achieved. 
Similar consideration hold for the experimental measure- 
ment of quantum Casimir forces [[[o| O, [l^] , but we focus 
here only on the classical regime. 

In this Letter we fully characterize the fluctuations of 
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FIG. 1: Two identical macroscopic plates normal to the x 
axis immersed in a fluctuating correlated medium. 



these thermally excited elastic forces. Generically, we 
consider in a space of dimension d, a fluctuating medium 
described by a scalar field <fi with an elastic energy density 
proportional to (V</>) 2 . For instance, 4> could represent 
the deviation of the director of a nematic liquid crys- 
tal pi from a uniform orientation (in the one-elastic- 
constant approximation) or the composition of a mix- 
ture at its critical point M. We consider two plates of 
lateral extension L separated by a distance H (see Fig. |l|) 
that impose 4> = on their boundary. We show that the 
fluctuating force is characterized by a Gaussian distri- 
bution for d > 1 centered on the usual Casimir force 
(F) ~ ksTL^ 1 1 H d . We characterize the corresponding 
dispersion by the "noise-over-signal ratio," i.e. the width 
AF of the distribution over the mean (F) . By calculating 
AF we show that this ratio scales as 



AF 



H 



(1) 



In the above result a corresponds to a microscopic length 
scale at which either the continuum description breaks 
down or the boundary conditions are not efficiently set 
by the plates Q . We also consider the dynamics of the 
field and the role of the finite temporal resolution r of 
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a measurement apparatus. If r a is the microscopic time- 
scale associated with the microscopic length a, we show 
that the noise to signal ratio is significantly reduced if the 
apparatus' resolution induces some additional averaging 
(i.e. r > r a ). More precisely, the average is still (F), 
but the noise to signal ratio now scales as: 



Re [<H 




(Fr) 



H 



(2) 



The relevance of the above results to experimental situ- 
ations is discussed at the end of the paper. 

We consider an elastic Hamiltonian of the form 



V4>(R) 



(3) 



with R = (x, r), and x the coordinate in the direction 
normal to the plates. We assume that the plates im- 
pose Dirichlet boundary conditions, i.e., <f>(x = 0,r) = 
(j)(x = H, r) = 0. The whole system is in contact with a 
heat bath imposing the temperature T. For this minimal 
model, the free energy !F(H) of the system can be exactly 
calculated as a function of the interplate distance H || . 
The Casimir force (F(H)) = —dF/dH, i.e. the ensemble 
average of the fluctuation-induced force, is given by the 
universal formula 



(F) = -k B T A d 



L 



d-l 



(4) 



where the prefactor A d = T(d/2)C,(d) / (Air) d / 2 depends 
only on the spatial dimension d. 

We shall pursue another route to calculate the Casimir 
interaction, which yields a clear intuitive picture of the 
origin of the fluctuation-induced force, and permits eval- 
uation of its fluctuations. We use the stress tensor asso- 
ciated with the field 0: T l0 = £5 l0 - (di<j>)[d£ /d{d 5 4>)}, 
where £ = \K{\/cf)) 2 is the energy density associated 
with Eq. (||). We first restrict our attention to the ef- 
fect of the interplate medium on the plate located at 
x = H . For a given conformation of the field, the pro- 
jection F < (H) on the x-axis of the elastic force exerted 
on the inner side of this plate is given by the integral of 
T XXl which reduces to 



F<(H) 



K 

y 



d*- X T[d x <f>{H,T)Y 



(5) 



since <f> vanishes all along the plate. Note that this 
quantity is always positive, therefore the field always 
pushes on the plate. Fourier transforming the field in 
the transverse direction: <f)q(x) = jd x v <j>(x, r) e _lq r 
yields H = £ q £q and F<(H) = £ q / q (JT) with 

£ q = ^L 1 - d Kfdx[<f^{x)(-B^+^)4>- q \ and / q (tf) = 
^L 1 ^ d K\d x (j)c l (H)\ 2 . These expressions suggest an in- 
teresting picture: each Eq can be understood as the 
energy of a three-dimensional string parameterized by 



FIG. 2: Typical configuration of a fluctuating string pinned 
at x — and x = H. The dashed cylinder sketches the 
effective cage imposed by the confining quadratic potential of 
curvature A"g 2 /L d_1 . 



[x, Re(0 q ),lm((/>q)]; each of these virtual strings has a 
line tension proportional to the rigidity K of the medium 
and is confined around the x-axis by a harmonic poten- 
tial of stiffness proportional to q 2 . These lines are pinned 
at x — and x — H as a consequence of the Dirichlet 
boundary conditions (see Fig. ||) and the force they ex- 
ert on the plate located at x = H is precisely fq(H). 
The total partition function of the medium between the 
plates is the product Z 
independent strings 



ri| q |<i/a Z q of those of the 



,0 q (ff)=O 

„ q , ^ q (x) e - B -W/ fcBT , (6) 

'0q(O)=O 

where a -1 is the high wavevector cutoff. The number 
of non-interacting strings describing our system is from 
simple mode counting: 

d-l 



N a 



L 



(7) 



Determining the distribution of the elementary forces 
fq{H) is straightforward. Since the free energy associ- 
ated with the string fields q is quadratic, the distribu- 
tion of dx4><i{H) is a Gaussian with zero mean and vari- 
ance d x d y Gq(H,H), where Gq(x,y) = (^ q (x)^)_ q (y)) is 
the two-point correlation function JL8| : 



G^x,y)=L a 



_j fceT sinh [q(H — u)] sinh^v) 



K 



q sm\i(qH) 



(8) 



where u = max(ir,y) and v = min(x,y). Since the dis- 
tribution of d x (f)q is Gaussian, the fluctuating elementary 
forces follow a ^-distribution: 



= f 



0(f) 



exp 



where 9 is the Heaviside step function, and 



fc R T 



q cothlqH) 

na 



(9) 



(10) 



The cutoff a" 1 appears in Eq. ( pTj[ ) because of the discon- 
tinuity of the derivative of G q (x, y) |l9| , as a manifesta- 
tion of the importance of the short wavelength fluctua- 
tions. From Eq. |^ follows that the elementary forces un- 
dergo large fluctuations: their standard deviation, given 
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by 



(11) 



compares with their average value. 

We can now consider the whole space consisting of 
three independent subsystems delimited by the plates. 
The instantaneous force F(H) experienced by the plate 
at x — H is the difference between two independent 
stochastic variables F < (H) — F < (oo), where F < (oo) ob- 
viously describes the contribution of the outer medium 
x > H. Within the string picture F(H) = £l[/q(-H") - 
/ q (oo)] thus appears as the sum of N a independent vari- 
ables. 

Mean Casimir force: From Eq. (|lO|), the mean 
force (F) experienced by each plate is given by (re- 
placing the discrete sum by an integral and disre- 
garding the cutoff as the integral converges): (F) — 
ifc B T(L/27r)' i - 1 /ff d - 1 q[ (? - qcoth(qH)}, which yields 
the result quoted in Eq. (f|) above. We thus recover the 
usual expression of the Casimir force, classically obtained 
by differentiation of the free energy. 

The cancellation of the non-universal cutoff contribu- 
tion in the former sum is easily understood in terms of the 
strings picture: since a string with label q has a correla- 
tion length ~ g -1 , the strings with q ^> H^ 1 do not feel 
the presence of the second plate and their average con- 
tributions on both sides of the plate cancel exactly. Con- 
versely, the force imbalance for each string with q < l/H 
scales as —k-sT/H, and there are (L/H) d ^ 1 such strings, 
the product yielding the scaling form of the Casimir force 
inEq. %. 

Distribution of the fluctuation-induced force: The vari- 
ance of the force, (AF) 2 = (F 2 ) — (F) 2 , is the sum of the 
2N a elementary contributions (A/ q ) 2 from both sides of 
the plate: 



(AF) 2 = (k B T) 2 L d - 



jd-1. 



responsible for the force fluctuations, the scaling behav- 
ior of Eq. (|l3|) is expected to hold for any macroscopic 
external object with smooth shape at the scale a. The 
force fluctuations also include a subdominant universal 
part A'F, which can be exactly calculated as 



(27T) 



d-1 



q coth (qH) 



2 



,(12; 



where the integral must be limited to |q| < a 1 to be de- 
fined. When H 3> a, the leading behavior is given by the 
limit H > 00, i.e., (AF) 2 ~ (fceT) 2 ^' 1 J 



q] 2 q d 2 dq, which yields 



AF 



(13) 



where again N a is the number of transverse Fourier 
modes defined in Eq. (^) above. Contrary to the Casimir 
(average) force, the variance of the fluctuation-induced 
force is not universal and is intrinsically related to the 
physics ruling the interaction between the elastic medium 
and the immersed plates (through the cutoff). Since the 
short wavelength fluctuations of the q strings are mainly 



A'F = (dA d ) 



1/2 



k B T 
H ' 



(14) 



This cutoff independent part of the dispersion originates 
from the Njj = (L/Hy 1 ^ 1 strings with q < ff _1 , as in 
the case of the mean Casimir force. 

The y/N2 factor in Eq. ([TJ]) can easily be understood 
from the central limit theorem applied to the extensive 
system of 2N a non-interacting strings (for d > f), gener- 
ating a force with mean and standard deviation of order 
k-eT/a. The central limit theorem holds since for all the 
strings A/ q (F")/AF — > as N a — > 00. As a consequence, 
the distribution Vf of the force experienced by the plate 
at x = H is a Gaussian: 



V F {F = f) = 



1 



a/2ttAF 



exp 



(f-(F)Y 

2(AF) 2 



(15) 



with its mean and variance being given by Eqs. (^J) and 
dl2|), respectively. Only the average force felt by the 
macroscopic objects is universal, and the dispersion of 
the force is essentially controlled by the plates-medium 
microscopic interactions. 

Estimation of the Casimir force dispersion: Assuming 
that the instantaneous fluctuation-induced force can be 
measured with a perfect precision, the noise over signal 
ratio is thus given by Eq. ([fl), obtained from Eqs. (j|) 
and ( fl3"| ) using N a — (L/a) . Let us estimate it in 
a typical experimental situation, e.g., a nematic liquid 
crystal with a director strongly anchored normal to the 
plates. In this case, we have d — 3 and a ~ 1 nm. If 
the typical size of the plates is L ~ 10 /im and their 
separation is H ~ 100 nm, then the amplitude of the 
fluctuations is a hundred times larger than the mean force 
(F) ~ I pN. In order to lower the noise over signal ratio 
to 1, the size of the plates must reach 1 mm! 

In all the previous discussions we have implicitly as- 
sumed d > 1. It is interesting to notice the singularity 
of the one dimensional (ID) case. The net force experi- 
enced by the plate can no longer be interpreted as the 
result of a large number of independent contributions. A 
ID elastic medium is indeed equivalent to the sole q = 
string. Consequently the distribution of forces F{H) and 
F(oo) are directly given by Eq. (|^) taking the limit q — > 0. 
However, even if the distribution is no longer Gaussian, 
the scaling of the noise over signal ratio in Eq. (1) still 
holds in the limit d — > 1 . 

Measurement of fluctuation-induced forces and tem- 
poral resolution: It has just been shown that the equi- 
librium distribution of fluctuation-induced forces may 
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be extremely broad. Any experimental measurement 
will however provide a filtered signal which is averaged 
over the temporal resolution t of the apparatus. In 
order to estimate to what extent the fluctuating na- 
ture of fluctuation-induced forces may be experimen- 
tally revealed, a dynamical description is required. A 
precise description of the motion of the plates is be- 
yond the scope of this tentative approach. In a sim- 
ple picture, the measurement apparatus provides a sig- 
nal F T (t) = f_ x(i— t')F(t')dt'. Its response function x 
has a typical decay time r and is causal and normalized 
(J dt' x(t') — !)■ The limit r — > would correspond to 
a perfect apparatus which provides a signal completely 
described by Eq. (|l5|). 

The short wavelength excitations of the strings domi- 
nate the fluctuations of the net force as it has been previ- 
ously shown. When measuring the force, the apparatus 
averages over N T ~ r/r a independent processes, in which 
T a is the microscopic correlation time of the fluctuations 
associated to the modes with wavevector ~ a -1 in the x 
direction. Consequently, we expect the force dispersion 
AF to be lowered by a factor of l/y/N T , which leads to 
the result reported above in Eq. (^). 

In order to check this analysis, we have studied the 
simplest case of a local and dissipative dynamics for the 
4>q fields described by the Langevin equations 

7 a t q (x, t) = K [d 2 x - q 2 ] ^(x, t) + £ q (z, t), (16) 

subject to </> q (0,i) = </> q (iJ, i) — 0. In the above equa- 
tion £ q is a Gaussian white noise with zero mean and 
correlations chosen so as to insure thermal equilibrium, 
i.e. (^(x,t)^(x',t')) =2^k B T(27T) d - 1 S d - 1 (q+q')S(x- 
x')S(t — t'). This set of N a equations correspond to 
the Rouse dynamical equations for the N a pinned elastic 
stings undergoing thermal fluctuations |20|. The linearity 
of the Langevin equations above [Eq. ( Jl 6[ ) ] ensures that 
the 4>q's are Gaussian fields, so that it is straightforward 
to determine the statistical properties of the measured 
force. 

The finite temporal resolution r does not modify the 
mean force which remains given by Eq. (4), (F T ) = (F). 
The dispersion of the measured forces can be expressed as 
the product of the ideal expressions calculated previously 
by an attenuation factor: 



AF T = AF Y(t/t ), 



(17) 



where for this dynamics r a is given by r a = 7a 2 / ' K . The 
asymptotic scaling behavior of T(r/r ) does not depend 
on the specific form of the response function \: we have 
T(0) = 1 and T(s) - l/*Js in the limit s > 1. This 
simple model of a diffusive dynamics thus clearly cor- 
roborates the qualitative analysis that lead to Eq. (j^) 
above. 

Let us determine the effect of the finite temporal 
resolution for the experimental example quoted above 



(L = 10 fim, H — 100 nm, a — lnm). The order of 
magnitude of the relaxation time r a in an experiment 
involving a nematic liquid crystals is around 0.1 /is at 
T = 300K. Using a measurement device characterized 
by r ~ 1 ms such as an optical tweezer [Q , the noise 
over signal ratio is lowered from 100 to 1. 

Generally, Eq. || quantifies to what extend a slow mea- 
surement device is best suited for experimental observa- 
tion of the universal part of fluctuation-induced forces. 
From the above analysis it also follows that replacing the 
fluctuation-induced force between small objects by the 
simple Casimir average when studying their interactions 
and collective behavior may not be justified if these ob- 
jects are fast movers. 

One of us (RG) would like to acknowledge ESPCI for 
hospitality during his visit, and support through the Jo- 
liot visiting chair. 
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